Universal fluctuations in heavy-ion collisions in the Fermi energy domain 
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We discuss the scaling laws of both the charged fragments multiplicity n fluctuations and the charge 
of the largest fragment Z max fluctuations for Xe + Sn collisions in the range of bombarding energies 
between 25 A ■ MeV and 50 A ■ MeV. We show at E tab > 32MeV/A the transition in the fluctuation 
regime of Z max which is compatible with the transition from the ordered to disordered phase of 
excited nuclear matter. The size (charge) of the largest fragment is closely related to the order 
parameter characterizing this process. 



Theoretical description of the fragment production in 
heavy-ion (HI) collisions depends on whether the equilib- 
rium has been reached before the system starts fragment- 
ing. Possibility of the critical behavior associated with 
the transition from the particle evaporation regime at 
low excitation energies to the explosion of the hot source 
at about 5-10 MeV/nucleon cannot be excluded. Un- 
fortunately, this exciting possibility is difficult to study 
because all standard models and methods of characteriz- 
ing different phases and transitions of the nuclear matter 
in HI collisions assume an equilibrium mechanism of the 
fragment production. In this work, we shall apply new 
methods of the theory of universal fluctuations of ob- 
servables in finite systems Q to examine what can be 
said in a model independent way about the fragmenta- 
tion mechanism and the phase changement in HI colli- 
sions in the Fermi energy domain. Our analysis, which 
is independent of the assumption of the equilibrium in 
the fragments production process, uses the data of the 
INDRA multidetector system for Xe + Sn collisions at 
25 MeV < Ei ab /A < 50 MeV g-||. 

Several features of finite systems are important if one 
wants to study either the criticality or the distance to the 
critical point Jl[. These are : 

(i) The A-scaling of the normalized probability distri- 
bution P <m > [wi] of the variable m for different 'system 
sizes' < m > : 



< m > A P. 



= $(*(A)) 



< A < 1 



Z(a) = {m — m*)/ < m > l 



(1) 



(2) 



where < m > and m* are the average and the most prob- 
able values of m respectively, and $(^(A)) is the positive 
defined scaling function which depends only on a single 
scaled variable Z(A)- If the scaling framework holds, the 
scaling relation (1) is valid independently of any phe- 
nomenological reasons for changing < m > . The scal- 
ing domain in (1) is defined by the asymptotic behaviour 
of P <m> [m] when m — > oo and < m >^ oo, but z/&\ 



has a finite value. The A— scaling analysis is very robust 
and can be studied even in small systems if the proba- 
bility distributions P<m>[w] are known with a sufficient 
precision. In small systems however, the value of scal- 
ing parameter A may differ slightly from its asymptotic 
value (lj ; 

(ii) The tail of the scaling function 3>(z(a)) and the 
anomalous dimension g. 

All these features are related to the properties of the 
scaling function which characterizes the finite system. 

If the infinite system experiences a second-order phase 
transition, and if m is the extensive (scalar) order pa- 
rameter then [Q : 

(i) At the critical point, the corresponding finite sys- 
tem exhibits the 'first-scaling law' (A = 1) and the tail 
of scaling function for large positive Z(a) : $(#(a)) ~ 
exp(— Z(A))' ^ s characterized by a large value of the ex- 
ponent v = 1/(1 — g) > 2. 

(ii) The finite system exhibits the 'second-scaling law' 
(A = 1/2) in the ordered phase, and the first-scaling law 
in the disordered phase. In both cases, the tail of the 
scaling function is gaussian (v — 2). Close to the critical 
point, one may find also the cross-over phenomenon from 
the first-scaling to the second-scaling law by the contin- 
uous A-scaling law, again with the gaussian tail of the 
scaling function. 

If the parameter m is not singular at the transi- 
tion, then its probability distribution follows the second- 
scaling law with the gaussian tail. 

There are two generic families of the fragment produc- 
tion scenarios for which the second-order phase transition 
has been identified. The family of aggregation scenarios 
contains both equilibrium models like the Fisher droplet 
model, the Ising model or the percolation model, and 
off-equilibrium models like the Smoluchowski model of 
gelation. In these models, the average size of the largest 
cluster < s max > is the order parameter [Q,|| and the 
exponent r of the power-law cluster-size distribution at 
the critical point (r > 2) : n(s) ~ s~ T , is related to the 
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anomalous dimension as: g = l/(r — 1). The second 
family includes fragmentation scenarios and contains 
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FIG. 1. Different characteristics of the largest fragment 
charge probability distributions P[Z max \ for central Xe + Sn 
collisions, (a) P[Z max \ in the variables of the first-scaling law 
for Eiab/A — 25 (asterisks), 32 (crosses), 39 (triangles), 45 (di- 
amonds) and 50 (circles) MeV. (b) P[Z ma x] in the variables 
of the second-scaling law for Ei a b/A — 25 and 32 MeV. (c) 
The normalized second cumulant moment IC2 (= fi^/fti 2 ) of 
P[Zmax] is plotted as a function of In < Z ma x > 2 (= hiK 2 ) 
together with the statistical error bars. The lines A = 1 and 
A = 1/2 are shown to guide the eyes. 

various hybrids of the fragmentation - inactivation - bi- 
nary model j^]. In this family of models, the average 
fragment multiplicity < n > is the order parameter of 
the critical shattering process and the exponent r at the 



critical point (r < 2) is related to the anomalous dimen- 
sion as : g = r — 1. The order parameters for these two 
different scenarios are not only relevant in the study of 
phase changes in HI collisions but they are also measur- 
able. In this work, we shall investigate the patterns of 
the largest fragment charge distribution P[Z max ] and the 
charged fragment multiplicity distribution P[n] for differ- 
ent centrality conditions and bombarding energies in the 
Fermi energy domain, using the methods of the theory of 
universal fluctuations of observables in finite systems [EJ . 

Fig. 1 shows the A-scaling features of P[Z max \— dis- 
tributions for central Xe + Sn collisions at 25 MeV < 
Eiab/A < 50 MeV. In the experiment a great effort 
was done to well identify in charge the different frag- 
ments produced and especially the heaviest ones [p|,^0[. 
For each collision energy about 20000 events are taken 
into account in the present analysis. These events are se- 
lected with the experimental centrality condition : com- 
plete events (i.e., more than 80% of the total charge and 
momentum is detected) and Qfi ow > 7r/3. The lat- 
ter quantity is a global observable defined as the angle 
between the beam direction and the main emission direc- 
tion of matter in each event, which is determined from 
the energy tensor. It has been shown for the reactions in 
the Fermi energy domain that events with small Qfi ow 
are dominated by binary dissipative collisions Jll] , p^ ,p|. 
For events with little or no memory of the entrance chan- 
nel, fi ow is isotropically distributed. The upper part of 
the figure (Fig. la) shows that P[Z max \— distributions 
for Ei a b/A — 39,45,50 MeV can be compressed into a 
single curve in the scaling variables of the first-scaling. 
The distributions for 25 and 32 MeV, which show strong 
deviations with respect to this scaling curve both near 
the maximum and in the tail for large zr&) , can be com- 
pressed into another single curve in the variables of the 
second-scaling, as shown in Fig. lb. It should be stressed 
that we do not optimize the value of A because of the ex- 
perimental (number of events) and theoretical (smallness 
of the system) limits, but we study whether the data is 
consistent with one of the two limits : A = 1/2 and 1, 
which have a particular significance in the scaling theory 
of phase-transitions. 

More straightforward global measure of scaling fea- 
tures is provided by the cumulant moments : «i =< 
m >,K2 =< m 2 > — < m > 2 ,k^ =< m 3 > —3 < 
to 2 >< m > +2 < m > 3 , etc. In case of the A-scaling, 



normalized cumulant moments ■ *^^ A ' — 
independent of the system size < m > |]j 

= r( A 



IC { q a> = Kj/(ki)« a , are 
Log of the 

normalized cumulant moment /C = x> of P[Z max ] 

is plotted in Fig. lc versus log of < Z max > 2 , i.e. 
versus log of (ki) 2 . The data for different < Z max >, 
i.e. for different bombarding energies, should lie on a 
straight line if the A-scaling holds. The slope of this 
line gives the value of A. It is seen that the higher 
energy branch (E lab /A = 39,45,50 MeV) follows the 
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line A = 1 (the solid line), in agreement with Fig. la. 
The point for E^/A = 25 MeV is clearly off this line. 
The point for Ei a b/A = 32 MeV aligns along the line 
A = 1/2 (the dashed line) passing through the point for 
Eiab/A — 25 MeV, but is also close to the line A = 1. 
In this 'near-crossing' case, one has to investigate higher 
order moments, viz the whole probability distribution. 
Fortunately, Figs, la, lb show clearly that the collisions 
at 32 MeV I A belong to the branch A = 1/2. 




FIG. 2. The multiplicity distributions of charged fragments 
for central Xe + Sn collisions are plotted in the variables of 
the second-scaling law. Notation of data at different collision 
energies is the same as in Fig. 1. 

The pattern of charged fragment multiplicity distri- 
butions P[n] does not show any significant evolution 
with the bombarding energy (Fig. 2), and the data 
is perfectly compressible in the scaling variables of the 
second-scaling, i.e. the multiplicity fluctuations are small 
(k q ~ (^i) 9 ^ 2 ) in the whole studied range of bombarding 
energies. The scaling features of experimental P[Z max ] — 
and P[n]— probability distributions in Figs. 1 and 2 are 
complementary and allow to affirm that the fragment 
production in central HI reactions in the Fermi energy 
domain follows the aggregation scenario and exhibits the 
transition at E^/A > 32 MeV between the two phases 
of excited nuclear matter with distinctly different pat- 



terns of Z n 



-fluctuations. 



Assuming that the change of the fragmentation regime 
is controlled by the source excitation energy, one may 
ask to which extent similar scaling features and transi- 
tions can be seen in more peripheral collisions. For that 
purpose, we have investigated a sample of events selected 
with the experimental condition : complete events and 
Qfiow > 7r/18. This sample contains mostly events of 
semi-central binary collisions which keep some memory 
of the entrance channel as manifested by the presence 
of two contributions in the velocity distribution of Z max 
corresponding to the quasi-projectile and quasi-target [^) . 
Fig. 3a shows that Z max — distributions for Ei a b/A = 45 
and 50 MeV can be compressed into a single curve in the 
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FIG. 3. Different characteristics of the largest fragment 
charge distributions P[Z max ] for Xe + Sn collisions, (a) 
P[Zmax] for semi-central collisions in the variables of the 
first-scaling law for Ei a b/A = 25 (asterisks), 32 (crosses), 39 
(triangles), 45 (full diamonds) and 50 (full circles) MeV . (b) 
P[Z m ax] for semi-central collisions at Ei a b/A = 25, 32 and 39 
MeV in the variables of the second-scaling law. (c) P[Z max ] 
for semi-central collisions at Ei a i/A = 45 and 50 MeV and 
central collisions at Ei a b/A = 39, 45 and 50 MeV in the 
variables of the first-scaling law. The notation for central 
collisions is the same as in Fig. 1. 

25, 32 and 39 MeV show significant deviations with re- 
spect to this scaling curve, both at the maximum and 
in the large-Z(A) tail. These three distributions are 
shown in Fig. 3b in the variables of the second-scaling. 
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P[Zmax]— distributions for Ei a b/A — 25 and 32 MeV can 
be perfectly compressed into a single curve. The data for 
Eiab/A = 39 MeV is close to this curve but, neverthe- 
less, shows some deviations for Z(i/2) — 1- This data, 
which seems to be intermediate between the scaling lim- 
its A = 1/2 and A = 1, could indicate a continuous 
change of A in the transition region. There are presently 
not enough data to verify this possibility. We note how- 
ever, that an optimal data compression for Ei^/A = 32 
and 39 MeV would indicate A ~ 0.6 for the latter one. 

The scaling pattern of P[Z max ] in semi-peripheral col- 
lisions follows a similar general evolution with the bom- 
barding energy (the excitation energy) as seen in central 
collisions, except that the branch A = 1 in semi-central 
collisions starts at higher energies than in central col- 
lisions. Fig. 3c presents all Z max — distributions which 
in semi-central and in central collisions belong to the re- 
spective first-scaling branches. It is non-trivial that these 
probability distributions, which correspond to different 
selection criteria, collapse approximately into the unique 
scaling curve. This example should encourage further in- 
vestigation of the feasibility of semi-central events in this 
kind of studies. One should also note that the scaling 
functions in the branch A — 1/2 are essentially identical 
in semi-central and central collisions (compare Figs, lb 
and 3b). 

The characteristic feature of critical behavior is the 
anomalous tail of the scaling function with v > 2. Using 
all P[Z max ] shown in Fig. 3c, we have fitted the scaling 
function $(ztW) for zn) > by : a exp(— b(zm — Zq) v ), 
where zq is the estimate of the most probable value of the 
distribution and a, b, v are the fitting parameters. We 
have found : v = 1.6 ± 0.4, what is incompatible with 
typical values (3.5 < v < 6) in the critical region for 
aggregation scenarios Q . The same procedure for central 
and semi-central collisions in the second-scaling phase 
yields : v = 1.8 ± 0.4. Hence, the existing data do not 
show the critical bahavior in the transition region from 
ordered to disordered phase. 

In conclusion, we have applied the theory of univer- 
sal fluctuations in the finite systems jjl to the symmet- 
ric HI reactions in the Fermi energy domain. This the- 
ory provides rigorous methods to characterize critical 
and off-critical behaviors both in equilibrium and off- 
equilibrium finite systems. This is an important novel 
aspect of the present analysis. Convincing scaling behav- 
ior of the Z max — and n— distributions have been found 
in central collisions. The present analysis shows that the 
fragment production in central HI collisions at around 
the Fermi energy is governed by the aggregation scenario 
with < Z max > as the order parameter. The changement 
of the regime of Z max — fluctuations from the second- 
scaling at low energies to the first-scaling at higher ener- 
gies, with the gaussian tail of the scaling function in both 
regimes, is compatible with the transition from the or- 
dered phase to the disordered phase. The scaling curves 



determined from central and semi-central collisions are 
the same but the bombarding energy corresponding to 
change in the scaling behavior is higher for more periph- 
eral collisions, in according with the expected stronger 
contribution of non-equilibrium effects in these collisions. 
The critical region, if exists in the nuclear fragmentation 
process, should be searched for in the narrow window 
of bombarding energies close to : Ei a b/A ~ 32 MeV. 
Its signature would be the anomalous tail v > 2 of 
P[Z m ax\, in case of the second-order phase transition, or 
the double-hump shape of P[Z max \— distribution for the 
first-order phase transition p| . Future studies of symmet- 
ric HI reactions, performed in small steps of bombarding 
energies, will hopefully allow to distinguish between the 
cross-over phenomenon and the other two scenarios which 
invoke existence of the phase-transition. 
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